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ABSTRACT 

Several recent papers claim the detection of a near infrared Extragalactic Background 
Light (EBL) intensity at 1.25 - 4 /im that exceeds the integrated light of galaxies by 
factors of > 3. When combined with a claimed optical detection of the EBL at 0.80 fj,m 
the EBL excess emission spectrum has a discontinuity at ~ l/zm. This discontinuity 
has given rise to an interpretation in terms of ultraviolet radiation emanating from the 
first generation of massive stars at redshifts of 7 - 20 (so called Population III stars). 
The interpretation of the NIR excess emission as being of extragalactic origin depends 
crucially on the model used in the subtraction of the Zodiacal Light, the dominant 
foreground contaminant. We estimate the Zodiacal Light at 0.80 ^m using on the one 
hand the measurement by Bernstein, Freedman & Madore (2002b), with corrections 
for some omitted effects of atmospheric scattering and calibration, and on the other 
hand the model of Kelsall et al. (1998). There is in neither case any evidence for a 
step in the EBL at ~ 1/im. We emphasize that in order to avoid systematic effects it is 
essential to use the same Zodiacal Light model (Kelsall et al. 1998) for both the NIR 
(1.25 - 4 ^m) and optical (0.80 /im) data. We emphasize, however, that our analysis 
does not allow a statement on the overall level of the NIR EBL. The contribution of 
the Diffuse Galactic Light to the "EBL excess" emission is estimated. It is found to 
be significant at 3 - 4 /^m and should be carefully evaluated in future measurements 
which aim at detecting an EBL signal at the level of ~10 nW m^^ sr^^, i.e. at the 
level of the integrated light of (known) galaxies. 
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1 INTRODUCTION 

The Extragalactic Background Light (EBL) in the ultravio- 
let, optical, and infrared wavelength regions can potentially 
provide a measure for the total luminous matter in the Uni- 
verse since the epoch of star and galaxy formation; for a 
review see Longair (2001). 

Detection of a near infrared Extragalactic Background 
Light (EBL) component, in excess of the integrated light 
from known galaxies, has recently been announced by sev- 
eral research groups using the COBE/DIRBE (for a review 
see Dwek, Arendt, & Krennrich 2005) and the Japanese 
IRTS instruments (Matsumoto et al. 2005). A "first detec- 
tion" of the optical EBL at 0.30, 0.55, and 0.80 fim has 
been claimed by Bernstein, Freedman & Madore (2002a,b) 
(hereafter referred to as BFM02a,b). Combining the claimed 
detections in these two wavelength bands, i.e. at 1.25 - 4 /im 
and 0.8 fim, there appears to be an indication of a discon- 
tinuity at ~ l/im, the EBL brightness rising by a factor of 
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~ 4 from the short-wavelength to the long-wavelength side 
of this step (see e.g. Matsumoto et al. 2005). The potential 
importance of this spectral discontinuity stems from its sug- 
gested interpretation as signature of the first generation of 
stars (s.c. Population III stars) with their rest-frame emis- 
sion at A > 0.912/^m redshifted by a factor of 1 -|- z ~ 8 — 20 
to the near infrared (Santos, Bromm, & Kamionkowski 2002; 
Salvaterra & Ferrara 2003; Cooray & Yoshida 2004; Madau 
& Silk 2005; Dwek, Arendt, & Krennrich 2005). 

The reality of the NIR EBL excess at 1 - 4 ^m and 
the 1 ^m spectral step have been contested both on obser- 
vational (see e.g. Aharonian et al. 2006; Dwek, Arendt, & 
Krennrich 2005) and on theoretical grounds (see e.g. Madau 
& Silk 2005). The present author (Mattila 2003) has argued 
that, because of insufficient or incorrect atmospheric scat- 
tering corrections and calibration, the claim by BFM02 for 
the detection of the EBL at 0.30, 0.55, and 0.80 /im was 
premature. 

In the present paper we pose the question whether the 
BFM02 measurement can be used to set a sufficiently low 
upper limit to the EBL at 0.80 fim which, in combination 
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with the claimed NIR EBL values at 1.25-4 /im, would jus- 
tify the claim for a step near ~ l/im in the EBL spectrum. 
Because of the large intensity of the Zodiacal Light in com- 
parison with the EBL its accurate subtraction is crucial for 
any statement on the EBL and its spectral shape. Therefore, 
we examine in detail the underlying ZL determinations and 
models used in the above mentioned works. We will also ex- 
amine the contribution of the Diffuse Galactic Light (DGL) 
to the observed 0.8 - 4 /jm sky brightness. 

2 ZODIACAL LIGHT ESTIMATION 

The basic formula used in the above mentioned papers for 
extracting the NIR or optical EBL is the following: 

IeBL = I tot — IzL — IlSL — I DGL (1) 

where Itot is the total sky brightness as measured from 
space, IzL the Zodiacal Light (ZL) as estimated using ded- 
icated observations or a model of the interplanetary dust 
cloud, IisL is the Integrated Starlight (ISL) from unresolved 
stars evaluated from star count models, and Idgl is the Dif- 
fuse Galactic Light (DGL). 

The first two terms are large, ~300 nW m~^ sr~^ at 
1-2 /Lim, as compared to the EBL which according to e.g. 
Matsumoto et al. (2005) and Cambresy et al. (2001) is ~50- 
70 nW m~^ sr~^ , but has conventionally been considered 
to be much smaller, ~10 nW m~^ sr~^, and mainly due to 
the integrated light of galaxies only. While the ISL term is 
smaller, it is still several times larger than the EBL signal 
due to known galaxies. Thus, the demands for the accuracy 
of the hot and Izl determinations are extremely high. 

2.1 ZL estimated from the observations of 

Bernstein, Freedman &: Madore (2002a, b) 

In the case of the claimed EBL detections at 0.30, 0.55, and 
0.80 /Km by BFM02 the Izl term was derived from ground 
based observations of the depths of Fraunhofer lines between 
0.39 and 0.51 /im. Large uncertainties are caused to ground 
based night sky photometry by the corrections needed for 
the atmospheric scattered light and extinction. Also, the 
method necessitates to calibrate the extended source re- 
sponses of two different telescopes, a ground based and a 
space borne one, to the same absolute scaling. 

In a re-discussion of the BFM02 ZL observations Mat- 
tila (2003) showed that corrections had to be applied be- 
cause of the following errors or omissions: (1) incorrect 
aerosol albedo; (2) omission of ground reflectance; (3) omis- 
sion of DGL as source of atmospheric scattering; (4) incor- 
rect aperture-correction factor. All these corrections were 
negative in sign and thus had the effect of decreasing the 
ZL value. Their combined effect to 7zl(0.80 /um) amounted 
to -9.4 to -13.1 10"'' erg/s cm^sr A, see Table 1 of Mattila 
(2003). This correction is ca. 10 times larger than the accu- 
racy of ±1 10"" erg/s cm^sr claimed by BFM02a for their 
ZL value and ca. 4-6 times as large as their claimed EBL 
value of 2.2 10~^ erg/s cm^ sr A at 0.80 /xm. Each one of the 
corrections (1) (4) has large uncertainties caused by the ill- 
defined atmospheric, ground refiectance and telescope PSF 
properties. We estimate that the resulting correction has an 
uncertainty of ~ ±50% or ~ ±5 10-" erg/s cm'^sr A. In an 



Erratum paper Bernstein, Freedman & Madore (2005) dis- 
cussed a number of errors and omissions in BFM02. While 
some of the errors noted by Mattila (2003) were addressed, 
they did not yet address the four major corrections (1) - (4) 
as discussed above. The original ZL estimate of BFM02b is 
given in line 2 of Table 1 and the (Mattila 2003) estimate, 
including these corrections (1) - (4), in line 4 of Table 1. 

Another point of concern is how BFM02 handled the 
combination of the individual systematic errors: this was 
done by assuming a flat probability distribution for each con- 
tributing source of error. This method produces combined 
systematic errors which are by factor 2-3 smaller than 
the conventionally determined errors obtained by adding 
in quadrature the individual systematic errors. One specific 
systematic error, important for the present purpose of esti- 
mating Izl at 0.80 /xm, is caused by the reddening of the ZL 
relative to the Solar colour. BFM02a estimate the reddening 
correction factor per 1000 A to be C(A) = 1.044 with a sys- 
tematic error of [—0.014, +0.006]. Extrapolating Izl from 
the reference wavelength of 0.465 /im to 0.80 pm the system- 
atic error caused by the uncertainty in the ZL reddening 
correction thus becomes [— 4.1%, -|-1.75%]. Combining this 
in quadrature with the other systematic errors the total sys- 
tematic uncertainty of /zi(0.8/im) becomes [—4.4%, -|-2.4%] 
or [-3.0, -1-1.6] 10"" erg/s cm^sr A (line 3 of Table 1). 

2.2 ZL estimated using the Kelsall et al. (1998) 

model 

The Zodiacal Light model which has been adopted in many 
of the NIR EBL studies, e.g. Matsumoto et al. (2005) and 
Cambresy et al. (2001), has been constructed by Kelsall et 
al. (1998). This model uses the full C70BE/DIRBE data set 
at 10 wavelengths between 1.25 and 240 /im and relies essen- 
tially on the seasonal variations of the ZL. The IPD emis- 
sion dominates the sky brightness at high galactic latitude 
areas and mid IR wavelengths, A « 10 — 60/im, a situa- 
tion which helps to constrain the model parameters at these 
wavelengths. 

In order to investigate whether the different approaches 
used for the ZL estimation at A « 1— 4/im (modelling) and at 
A = 0.80;um (ground based observations, BFM02b) could be 
the reason for the ~1/Lim EBL spectrum step we have applied 
the Kelsall et al. (1998) ZL model to the Itot observation 
of BFM02a at 0.80/im. To do this we have extracted from 
the DIRBE Sky and Zodi Atlas (DSZA) the ZL value at 
J band (1.25/Lim) for the specific sky position (A - Aq, /3) 
and day of the year of the HST measurement of Itot by 
BFM02a. This ZL value, in units i)f M.Iysr^^, is given in 
line 2 of Table 2. The systematic error estimate given in 
brackets corresponds to the error of 15 nW/m^ sr"^ as given 
in Table 7 of Kelsall et al. (1998). To transform the tabulated 
DSZA value into a monochromatic intensity at 1.25 /xm a 
small colour correction is applied assuming for the ZL a solar 
temperature blackbody spectrum. 

To estimate the ZL intensity at 0.80 /um we adopt the 
Solar SED as given by Colina, Bohlin, R.C.& Castclh (1996), 
convolve it with the HST WFPC2 filter F814W transmis- 
sion curve using the SYNPHOT/calcphot procedure^, and 

^ http: / /www.stsci.edu/resources/softwcire_haTdwcire/stsdas/synphot 
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Table 1. Estimation of the Zodiacal Light and Extragalactic Background Light for the BFM02 Field at 0.8 //m using the ground-based ZL 
observations of Bernstein, Preedman & Madore (2002b) (BFM02), with corrections according to Mattila (2003) (M03). Systematic errors are 

given in lirackots. 





Quantity 


Unit 








Value 


Comment 


1 


74,(HST), 8000 A 


10-9 


erg/s 


cm^ 


sr A 


72.4 ±0.14, [±1.0] 


BFM02, with "1-ct" systematic error 


2 
3 
4 


/|^(0.80 ^m) 
(0.80 Mm) 
/|^(0.80/im) 


10-9 
10-9 
10-9 


erg/s 
erg/s 
erg/s 


cm^ 
cm^ 
cm^ 


sr A 
sr A 
sr A 


69.4 ±0.4,[-0.9, +0.8] 
69.4 ±0.4,[-3.0, +1.6] 
56.3-60.0 


BFM02, with "l-ir" systematic error 

BFM02, but with conservative systematic error 

BFM02 as corrected by M03 


5 
6 


-fsBi+DGiCO-SO/im) 
-ffiSi+DGiCO-SOMm) 


10-9 
10-9 


erg/s 
erg/s 


cm^ 
cm^ 


sr A 
sr A 


3.0 ± 0.4, [±0.9] 

3.0 ±0.8, [-1.9, +3.2] 


BFM02, with l-cr statistical and 
"l-cr" systematic errors 
BFM02, but with 2-ct statistical 



and conservative systematic errors 

7 7^g^_^QQ^(0.80/im) 10-9 erg/s cm^ sr A < (12.4 - 16.1) BFM02 with corrected ZL from M03 

8 ^■^EBL+DGLi^-^^l^'^) nW/m^ sr < (99 - 129) BFM02 with corrected ZL from M03 



Table 2. Estimation of the Zodiacal Light and Extragalactic Background Light for the BFM02 Field at 0.8/im using the Kelsall et al. (1998) 
(K98) Zodiacal Light model. Systematic errors axe given in brackets. 



Quantity 


Unit Value 


Comment 


1 /t^„t(HST), 8000 A 


10-9 erg/s cm^sr A 72.4 ± 0.14, [±1.0] 


BFM02, with "I-ct" systematic error 



2 


/^j^ (model), J band 


MJy 


sr-l 






0.151, [±0.0062] 


K98 model value from DSZA 


3 


Colour correction factor K 










0.985 


Assuming 6000 K blackbody 


4 


/-^ (1.25 Mm) 


MJy 


sr— 






0.153, [±0.0063] 


Monochromatic intensity 


5 


7^^(1.25 Mm) 


10-9 


crg/f 


i cm^ 


sr A 


29.4, [±1.15] 


Monochromatic intensity 


6 


(1.25 Mm) 


erg/s 


cm 


A 




44.6 


Monochromatic Solar flux 


7 


(0.80 Mm) 


erg/s 


cm^ 


A 




112.3 


Solar flux through WFPC f814 


8 


(0.80 Mm) /F^ (1.25 Mm) 
/|^(0.80Mm)//|^(1.25Mm) 










2.52 


Solar SED 


9 










2.26[±0.05] 


Reddened Solar SED 


10 


4^(0.80 Mm) 


10-9 


erg/s 


i cm^ 


sr A 


66.4[±2.6], [±1.5] 


ZL resulting from the K98 model 


11 


^BBL+BGL(0-80A'm) 


10-9 


erg/s 


i cm^ 


sr A 


6.0 ±0.1, [±3.2] 


I^„t(BSr£) minus /zi(K98) 


12 


^•-fBSi+i3GL(0-80Mm) 


nW/m^ sr 






48.0 ± 25.6 


A • /^^(HST) minus A • /zi,(K98) 



account for the ZL reddening relative to the Solar SED by 
applying a correction factor according to Eq.(22) of Lcincrt 
et al. (1998). The reddening correction 1.25 fim —> 0.80 fim 
of 1.116 is estimated to have a systematic error of ±0.023, 
i.e. ±20%. The resulting ZL brightness of 66.4[±2.6], [±1.5] 
at 0.80 nm is given in line 10 of Table 2. The first error in 
brackets corresponds to the Kelsall ct al. (1998) model un- 
certainties, the second one to the 20% error adopted for the 
ZL reddening correction. 



2.3 The residual signal, EBL -|- DGL, after 
subtraction of the ZL 

BFM02a give for their HST measurement of Itot at 0.80 
yum the value 72.2±0.14 [±1.0] 10-9 erg/s cm^srA(line 1 

of Table 1 and Table 2). Their method of treating system- 
atic errors leads to the paradoxical result (line 5 of Ta- 
ble 1) that subtracting the value of /z_L(0.8Mm) = 69.4 
±0.4[-0.9, +0.8] 10-9 erg/s cm^srA leaves as difference 
the value Iebl+dgl (O.Sfim) = 3.0 ±0.4[±0.9] 10-9 erg/s 
cm^sr A , i.e. the systematic error has decreased from its level 
for /tot(0.8/im) alone. We give in line 6 of Table 1 a more 
conservative error estimate with 2a statistical and quadrati- 
cally added systematic errors. The systematic errors of ~ ±5 



10-9 erg/s cm^sr A caused by the four corrections discussed 
by Mattila (2003) and in Sect. 2.1 above have not boon in- 
cluded here. Instead, lines 7 and 8 give the estimated upper 
limits of Iebl+dgl and A • Iebl+dgl, reflecting the con- 
clusion of Mattila (2003) that the claim for a detection of 
the EBL by BFM02 viras premature. 

Subtraction of the Kelsall et al. (1998) Izl esti- 
mate from Itot results in the values of Iebl+dgl and 
^ • Iebl+dgl as given in lines 11 and 12 of Table 2. The 
systematic errors of hot and Izl have been added in quadra- 
ture. 



Fig. 1 shows the A • Iebl+dgl values and upper limits. 
They include the IRTS Near-Infrared Spectrometer (NIRS) 
values of Matsumoto et al. (2005) between 1.4-4 Mm and 
the COBE/DIRBF, values at 1.25, 2.2, 3.5, and 4.9 Mm as 
compiled in Dwek, Arendt, & Krennrich (2005) and Hauser 
& Dwek (2001). Our two estimates at 0.8 Mm are shown as 
an upper limit (box), based on the BFM02b ZL estimate as 
revised by Mattila (2003), and as an open circle with error 
bars, based on the ZL estimate from the Kelsall et al. (1998) 
model. 
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Table 3. Estimation of the diflfuse Galactic surface brightness due to scattered starlight, lDGL(sca.), and non-equilibrium emission or fluores- 
cence by dust, /i3<3x,(em). 



Field 
(1) 


l,b 
deg 

(2) 




JV(H) 
1020 cm-2 

(3) 


Av 

mag 

(4) 


A 

fini 
(5) 


nW/m sr 
(6) 


Ax 

mag 
(7) 


(8) 


A/^Q^(sca) A/^Q^(em) 

nW/m sr 
(9) (10) 


^^DGL (sca-|-em) 
(11) 


BFM02 


206.6, 


-59.8 


1.5 


0.080 


0.55 
0.80 


660" 

867" (900) 


0.080 
0.048 


53 
42 


21^' 

17 (8.0'=) 




21 
17 


Matsumoto 
et al. (2005) 


80-90, 


42-48 


1.3 


0.069 


1.25 
1.65 
2.2 


1060*^ 
1110'' 
655^- 


0.019 
0.012 
0.0075 


20 
13 
4.9 


7.9 
5.2 
1.9 


0.3 
1.2 
3.8 


8.2 
6.4 
5.7 



3.5 5.0 



" Mattila 1980 

^ derived from I DGL{sca)/ 1 (lOO^im) in L 1642 (Laurcijs, Mattila & Schnur 1987) 
effective value to be applied to Iebl+dgl oi BFM02 
Lehtinen & Mattila (1996) 



3 THE CONTRIBUTION BY THE DIFFUSE 
GALACTIC LIGHT 

We have so far presented in Tables 1 and 2 and Fig. 1 the 
combined contribution of the EBL and the DGL. Since Mat- 
sumoto et al. (2005) do not make any statement on the DGL 
they appear to have included it into their nominal "EBL" 
values. For the COBE/DIRBE based values at 1.25, 2.2, 
3.5, and 4.9 shown in Fig. 1 the DGL has been considered 
to different degrees. The following discussion in this section 
is concerned with the DGL contribution to the Matsumoto 
et al. (2005) 1.4 - 4/im and the BFM O.SOjum values. 

BFM02a used a general model of the ISRF and dust 
distribution to estimate the DGL (scattered starlight) con- 
tribution in their target field at 0.80 /xm: Idgl = 0.8 10~® 
erg/s cm^srA. Mattila (2003) revisited the model estimate 
and argued that the dust column density toward the target 
field had been underestimated by a factor of ~ 3 and esti- 
mated the DGL intensity at 0.8/im to be Idgl = 2.3 10~^ 
erg/s cm^sr A corresponding to 18 nW/m^sr. However, after 
including a correction caused by the fact that the DGL has a 
Fraunhofer line spectrum resembling the ZL spectrum Mat- 
tila (2003) found that an effective Idgl correction of ~ 1 
10"^ erg/s cm"^sr Aor 8 nW/m^sr shall be subtracted from 
the observed Iebl+dgl value of BFM02. 

At optical and NIR wavelengths up to 1.5 /xm the dif- 
fuse Galactic surface brightness is primarily due to Galac- 
tic starlight scattering off interstellar dust grains. At longer 
wavelengths (~ 2 — 5/im), because of decreasing optical 
depth of dust {A\ < 0.15 Ay), the scattered light contri- 
bution rapidly decreases. However, recent COBE/T)IRBE 
and Spitzer observations have shown that there is at NIR 
wavelengths, A > 2/im, a substantial diffuse Galactic sur- 
face brightness component (Bernard et al. 1994; Arendt et 
al. 1998; Arendt & Dwek 2003; Flagey et al. 2006) associ- 
ated with the widely distributed low density dust medium. 
Such a NIR continuum emission had been previously de- 
tected in reflection nebulae by Sellgren, Werner & Diner- 
stein (1983). It can be due e.g. to non-equilibrium emission 
by very small transiently heated grains or to fluorescence 



emissions by PAHs or very small grains; for a discussion see 
Flagey et al. (2006). 

We will estimate the diffuse Galactic emission at the 
different wavelengths and for the two regions of sky tar- 
geted in BFM02 and Matsumoto et al. (2005) observations. 
The line-of-sight dust column density can be estimated us- 
ing the 21-cm hydrogen column density map of Kalberla et 
al. (2005). For this map the best state-of-the-art stray radi- 
ation corrections have been applied, important for the high 
latitude flelds with weak HI emission. The HI column den- 
sities for the two fields are 1.5 and 1.3 lO^o cm"^ as given in 
column(3) of Table 3. Using the standard N(U)/E{B - V) 
ratio (Bohlin, Savage, & Drake 1978) and Rv = 3.1 they 
correspond to the line-of-sight extinctions of Ay = 0.080 
and 0.069 mag, respectively. We assume that the scatter- 
ing properties of the dust grains, i.e. albedo a and forward 
scattering parameter g, are approximately constant between 
0.55 and 2.2 /xm (for justification see Lehtinen & Mattila 
1996). Then, the scattered light intensity at different wave- 
lengths is, for an optically thin line of sight towards high 
galactic latitudes, in a first approximation proportional to 
the product /(ISRF) x Ax- We give in columns (6)-(8) of 
Table 3 our estimates for /(ISRF) and Ax- The wavelength 
dependence of Ax corresponds to Rv = 3.1 (Mathis 1990). 
For /(ISRF) we have adopted at 1.25 — 2. 2/im the mean sky 
brightness values as obtained from the COBE/DIRBE Zodi 
Subtracted Mission Averaged (ZSMA) maps; see Lehtinen 
& Mattila (1996). 

The scattered light intensity, /DGL(sca), has been mea- 
sured in the translucent high latitude (6 = —36.7 deg) 
cloud L 1642 as a function of optical thickness for five 
wavelengths between A = 0.35 - 0.55 /im (Laureijs, Mat- 
tila & Schnur 1987; Mattila 1990). At smaU optical depths 
we obtain from these observations at 0.55 /im the relation 
/i3Gi(sca)//(100^m) = 3.0 10"^ erg/s cm^srA/MJy sr'^ 
which for /(100/im)/iV(HI) = 0.85 MJy sr'VlO^" cm'^ 
(Boulanger et al. 1996) corresponds to /_DGL(sca)/N(i//) = 
2.55 10"* erg/s cm^sr A /lO^" cm'^. Thus, for the BFM02 
fleld we obtain at 0.55 /xm the estimate A/^£;^(sca) = 21 
nW/m sr as given in column (9) of Table 3. The A/^Gi(sca) 
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1 I — ' — ' — ^ ' ' ' 1 

1 

Wavelength [/im] 

Figure 1. Observed values and upper limits for the Extragalactic 
Background Light in the red and near infrared A = 0.8-5 /im. The 
two estimates at 0.8 jxm. derived in this paper are indicated by an 
open circle and error bars (ZL model of Kelsall et al. 1998) and a 
rectangular box (BFM02 ZL observations with corrections accord- 
ing to Mattila 2003). The open circles with error bars at 1.4-4 /im 
are the IRTS Near-Infrared Spectrometer (NIRS) values of Mat- 
sumoto et al. (2005). The COBB/DIRBE values at 1.25, 2.2, 3.5, 
and 4.9 //m are indicated by big solid symbols: squares (Dwek 
& Arendt 1998; Arendt & Dwek 2003), triangles (Cambresy et 
al. 2001), octagons (Wright 2001), and encircled octagons (Gor- 
jian, Wright & Charry 2000). The integrated light from known 
galcixies is shown by small solid symbols: hexagons (Madau & 
Pozzetti 2000), triangles (Totani et al. 2001) and squares (Fazio 
et al. 2004). The estimated lower and upper limit for the total 
light of galaxies is shown by bars at 0.61, 0.81, 1.25, and 2.2 ^m 
(Totani et al. 2001). Our estimate for the Diffuse Galactic Light, 
^-DGi (sca-t-em), is shown as a dashed line according to column 
11 of Table 3 



values for the other wavelengths in Table 3 are then derived 

by scaling this A/^Qj^(sca) at 0.55 /im with the XIjsrf >^ 
values given in column (8). The value thus obtained for 
the BFM02 field at 0.80 /im is 17 nW/m^sr which very 
closely corresponds to the model estimate of BFM02a as 
corrected by Mattila (2003). These optical DGL estimates 
are in general terms also confirmed by the surface photome- 
try of L 1780 (Mattila 1979) which stretched over the larger 
wavelength range of 0.35 - 0.755 /im but lacked the smallest 
optical depth range important for the present purpose. It 
can be seen from Table 3 column (9) that XI^qj^{scsl) drops 
rapidly at A > 1.5/im. 

We use the results of Flagey et al. (2006) to estimate 
the XIdgl values at 1.25 - 4.9 /im for the Matsumoto et al. 
(2005) field. In their Fig. 9 they give the diffuse Galactic 



brightness at 2 - 14 fim for A''(H) = 10^^ cm~^. From this 
figure we read the loGLitot) — looLisca) + /oci (emission) 
values for 2.2, 3.5, and 4.9 /xm and they are given in column 
(11) of table 3. In order to estimate /ijoz, (emission) at 1.25 
and 1.65 /im we use the Flagey et al. (2006) model for the 
NIR continuum which they represent with a grey body of 
colour temperature 1100±300 K. The values thus obtained 
for A/oGi (emission) are given in column (10) of Table 3. 
At these two wavelengths Jdgl (sea) already dominates over 
/dgl (emission). It can be seen from Table 3 that the total 
XIdgl in units of nW/m^sr varies by a factor of ~2 over the 
NIR wavelength range, A = 1.25 - 4.9 jjxn. At 0.80 nm it is 
~ 2 times as large as at 1.25 //m. However, after including 
the correction caused by the similarity of the DGL and ZL 
Praunhofer line spectra its effective value, to be subtracted 
from the BFM02 XI^gj^_^_jjGL value at 0.80 /im, becomes 8 
nW/m^sr, i.e. very similar to the value at 1.25 /tm. 

The DGL contributions are, in general, substantially 
smaller than the XI^bl+dgl values of BFM02 and Mat- 
sumoto et al. (2005) shown in Fig. 1. The values are, how- 
ever, closely similar to the integrated intensities derived 
from galaxy counts. Therefore, if the true NIR EBL sur- 
face brightness turns out to be of this order of magnitude 
the DGL contribution cannot be neglected. Rather, it will 
be the ultimate obstacle, even for space borne measurements 
from outside the Zodiacal cloud, which must be accurately 
known and subtracted before arriving at a credible measure- 
ment of the EBL. 



4 DISCUSSION 

4.1 Is there an ~1 /im step in the EBL spectrum? 

The interpretation of the NIR EBL excess as signature of 
Population III stars hinges on the reality of the claimed 
large discontinuity between 0.8 and 1.2-1.4 /im. The values 
(upper limits) derived in Section 2 and shown in Fig. 1 for 
^Iebl+dgl at 0.8 /tm are significantly larger than the origi- 
nal BFM02a value of 3.0 10"® erg/s cm^ sr A or 24 nW/m^sr 
on which this conjecture was based. 

The large range of the estimates at O.SO/tm reflects 
directly the differences and uncertainties of the values of 
/zl(0.80 /im) as derived from the Kelsall et al. (1998) model 
on the one, and from the BFM02b ground based observa- 
tions on the other side. In both cases the ZL value has been 
extrapolated from another wavelength, 1.25 /tm in the for- 
mer and 0.465 /tm in the latter case, using a reddened Solar 
SED to transform from Solar to ZL colour. 

When assessing the reality of the ~1 /tm step in the EBL 
spectrum it is obvious that the subtraction of the foreground 
ZL has to be done using consistent ZL values, both on the 
short and on the long wavelength side of the step. With this 
in mind we have applied the Kelsall et al. (1998) model and it 
can be seen from Fig. 1 that there is, within the error limits, 
no evidence for the step. Errors in the basic assumptions or 
parameter values of the Kelsall et al. (1998) model would 
influence the values at 0.80 /tm and at > 1/tm in a similar 
way and are not expected to artificially create or destroy a 
large spectral step, if existent, between these wavelengths. 

When applying the BFM02b measurement of Izl an 
extrapolation of their measured ZL value from 0.465 /im to 
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0.80 fim is needed. This introduces an uncertainty of similar 
size as the extrapolation of the Kelsall et al. (1998) model 
from 0.80 to 1.25 fj,m.. Given the corrections needed for the 
BFM02b ZL measurement and adopting conservative sys- 
tematic error estimates as discussed in Sect. 2.1 above, there 
is no evidence for a downward step in the EBL spectrum 
from 1.25 to 0.80 fj.m. 

The Diffuse Galactic Light correction as discussed in 
Sect. 3 is of similar strength both below and above 1 /jm 
and it does not, therefore, contribute to the possible discon- 
tinuity at ~1 /jm. 

The author wishes to emphasise that he does not advo- 
cate the (very) large Iebl+dgl value at 0.80/im, apparently 
resulting from the discussion in this paper, as anything else 
but an upper limit. Rather, also the value derived using the 
Kelsall et al. (1998) ZL model, which appears to be different 
from zero, should be considered only as an upper limit. 

4.2 Remarks on the near infrared EBL 

measurements 

The analysis presented in this paper does not allow conclu- 
sions concerning the overall level of the NIR EBL at 1 - 4 /xm. 
Given the systematic uncertainties of the ZL subtraction a 
direct photometric detection of the NIR EBL appears very 
difficult. 

The Kelsall et al. (1998) model appears to give a good 
representation of the ZL thermal emission at mid and far 
IR wavelengths where the ZL dominates the sky brightness. 
However, at the near IR wavelengths, l/um < A < 5/xm, the 
Zodiacal component is weaker and its modelling relies on the 
mid IR results. As stated by Kelsall et al. (1998) their ZL 
model is not unique but the predicted values depend on the 
assumed models for the IPD cloud geometry and dust pa- 
rameters. FYirthermore, the NIR ZL is caused by scattering 
instead of thermal emission which means that dust parame- 
ters different from those at the mid IR wavelengths have to 
be used. 

An alternative ZL model has been presented by Wright 
(1998) and has been applied in the Wright & Reese (2000), 
Wright (2001), and Gorjian, Wright & Gharry (2000) anal- 
yses of the DIRBE data. The resulting EBL values are by 
~33, 8, and 4 nW/m^sr smaller at the J, K, and L bands, 
respectively, as compared to the results based on the Kelsall 
et al. (1998) model. This does not, however, substantially 
change the dilemma with an unexpectedly large NIR EBL 
value. 

Dwek, Arendt, & Krennrich (2005) have presented ar- 
guments suggesting that the Matsumoto et al. (2005) EBL 
excess emission at A = 1.4 — 4 /im may be due to insufficient 
ZL subtraction. This possibility is especially suggested by 
the practically identical SEDs of the ZL and the "EBL ex- 
cess" over the whole wavelength range of the observations, 
see Fig. 6 of Dwek, Arendt, & Krennrich (2005). 

The "EBL excess" emission as derived from the IRTS 
data by Matsumoto et al. (2005) contains the contribution 
by the Diffuse Galactic emission. The DGL has been con- 
sidered in several of the analyses of the COBE/DIRSE data 
(sec e.g. Arendt ct al. 1998; Arendt & Dwek 2003) but its 
separation from the much larger component due to the Inte- 
grated Starlight has been problematic because of the large 
field of view of DIRBE. The recent Spitzer observations 



(Flagey et al. 2006) have offered a better estimate for the 
2-5 /im DGL and they suggest that it is substantially larger 
than the values previously found by Arendt et al. (1998). Al- 
though the DGL level remains insignificant over the wave- 
length range 1.25-2.5 /im this is no longer the case at 3.5 
ytim. At this wavelength the sum of the integrated light from 
resolved galaxies (Fazio et al. 2004) and the DGL (Table 3) 
amounts to ~10.4 nW/m^sr which is within its estimated 
error limits of ^ ±2 nW/m^sr fully compatible with the 
Matsumoto et al. (2005) EBL value of 14.4 ± 3 nW/m^sr. 
In addition, similarly cis at the shorter wavelengths (Totani 
et al. 2001), some contribution from the unresolved galaxies 
should be added (Kashlinsky ct al. 2005; Savage & Oliver 
2005) bringing the two estimates even closer to each other. 

Recently, an upper limit to the EBL flux around 1 - 2fim 
has been announced by Aharonian et al. (2006) using HESS 
observations of intergalactic absorption of 7-ray emission of 
blazars: /ebl < (14±4) nW/m^sr. Thus, these obser- 
vations would suggest that more than two thirds of the EBL 
in the NIR band is resolved into individually detected galax- 
ies. The HESS result depends, however, on the adopted TeV 
spectra of the blazars. Therefore, direct photometric obser- 
vations of the optical/NIR EBL are still urgently needed. 



5 CONCLUSIONS 

We have scrutinised the observational evidence for the 
claimed discontinuity at ^1/jm in the spectrum of the 
Extragalactic Background Light. We have also estimated 
the contribution of the Diffuse Galactic emission to the 
"EBL excess" . Our conclusions are the following: 

(1) The reality of the ^Ifiui step hinges on the claimed 
detection of the EBL by BFM02 at 0.80 ixm. We find that 
after applying corrections to the analysis of their Zodiacal 
Light measurement only an upper limit to the EBL can be 
set which does not imply any spectral step between 0.80 
and 1.25 /im. 

(2) An alternative estimate of the Zodiacal Light at 0.80 
/xm is obtained using the Kelsall et al. (1998) model. Again, 
the upper limit set to the EBL at 0.80 pm does not warrant 
any spectral step between 0.8 and 1.25 /jm. 

(3) An estimate for the contribution by Diffuse Galactic 
emission is obtained from recent Spitzer measurements by 
Flagey ct al. (2006). It is found that the sum of the Diffuse 
Galactic emission and the integrated light of resolved 
galaxies can explain the claimed EBL signal in the 3.5 yum 
window where the ZL contamination is at minimum. 
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